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Abstract 

The finite temperature Casimir free energy, entropy, and internal energy are considered anew for 
a conventional parallel-plate configuration, in the light of current discussions in the literature. In 
the case of an "ideal" metal, characterized by a refractive index equal to infinity for all frequen- 
cies, we recover, via a somewhat unconventional method, conventional results for the temperature 
dependence, meaning that the zero-frequency transverse electric mode contributes the same as the 
transverse magnetic mode. For a real metal, however, approximately obeying the Drude dispersive 
model at low frequencies, we find that the zero- frequency transverse electric mode does not con- 
tribute at all. This would appear to lead to an observable temperature dependence and a violation 
of the third law of thermodynamics. It had been suggested that the source of the difficulty was 
the behaviour of the refiection coefficient for perpendicular polarization but we show that this is 
not the case. By introducing a simplified model for the Casimir interaction, consisting of two 
harmonic oscillators interacting via a third one, we illustrate the behavior of the transverse electric 
field. Numerical results are presented based on the refractive index for gold. A linear temperature 
correction to the Casimir force between parallel plates is indeed found which should be observable 
in room-temperature experiments, but this does not entail any thermodynamic inconsistency. 

PACS numbers: ll.lO.Wx, 05.30.-d, 73.61.At, 77.22. Ch 



* Electronic address: |johan.hoye@phys.ntnu.no| 
Electronic address: liver. h. brevi kQmtf .ntnu.nol 
^Electronic address: Ij an . b . aar set h@ mt f . nt nu . no| 
^Electronic address: niilton@nhn.ou.edui 



I. INTRODUCTION 



In spite of the numerous treatises on the Casimir effect during the past decade — for 
some books and review papers see, for instance, Milton Mostepanenko and Trunov 0], 
Milonni Plunien et ah fl, Bordag et ah |5]— it is somewhat surprising that such a basic 
issue as the temperature dependence of this effect is still unclear and has recently given rise 
to a lively discussion. This issue is not restricted to the case of curvilinear geometry, but is 
present even in the simplest conventional geometry of two parallel metal plates separated by 
a gap of width a. Thus Klimchitskaya and Mostepanenko in their detailed investigation p, 
and also Bordag et al. 01, and Fischbach et al. j^, have argued that the Drude dispersion 
relation for a frequency-dispersive medium leads to inconsistencies in the sense that the 
reflection coefficient r2 for perpendicular polarization (the TE mode) becomes discontinuous 
as the imaginary frequency C = —iuj goes to zero. As is well known, the Drude dispersion 
relation reads for imaginary frequencies 

where ujp is the plasma frequency and u the relaxation frequency. (Usually, u is taken to be 
a constant, equal to its room-temperature value.) The mentioned authors, instead of the 
Drude relation, give preference to the plasma dispersion relation, since no such discontinuity 
is then encountered. (In Ref. the plasma relation together with the so-called surface 
impedance approach is argued to be the method best suited to describe the thermal Casimir 
force between real metals.) The plasma relation is 



si^C) = l + ^. (1.2) 

The arguments in Refs. 0, 0, S, 1^ are interesting, since they raise doubts not only about 
the applicability of the Drude model as such, but even more, doubt about the applicability 
of the fundamental Lifshitz formula at low temperatures (see, for instance, Ref. |lO|). 

The essence of the problem appears to be the following: For a metal, does the transverse 
electric (TE) mode contribute to the Casimir effect in the limit of zero frequency, corre- 
sponding to Matsubara integer m = 0? It is precisely for this mode that the purported 
discontinuity of the reflection coefficient r2, mentioned above, can occur. The problem is 
most acute in the high T regime (the m = contribution becomes increasingly important 
as T increases), but is present at moderate and low temperatures as well. The conventional 
recipe for handling the two-limit problem for a metal, n = y/e — > cxd, m — > 0, has been to 
take the limits in the following order: 

1. Set first e = oo; 

2. then take the limit m = 0. 

This way of proceeding was advocated in the early paper of Schwinger, DeRaad, and Milton 
(we will call it the SDM prescription), and was followed also in one of the recent papers 
by some of the current authors 12], and in Milton's recent book It seems to escaped 
recent notice that the physical basis for this prescription, namely the necessity of enforcing 



the correct electrostatic boundary conditions, was explicitly stated in Ref. [11 
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Bostrom and Sernelius 3| seem to have been the first to inquire whether this prescription 
is right: They argued that in view of a realistic dispersion relation at low frequencies the 
m = TE mode should not contribute. And three of the present authors arrived recently at 
the same conclusion, in two papers dealing with the case of two concentric spherical surfaces 

aiii. 



The Bostrom- Sernehus paper gave rise to a heated debate in the literature 0,0,113,113 
on the role of the m = TE mode for a metal. The advent of accurate experiments in recent 
years, by Lamoreaux 0, Mohideen et al. JE, H, 2^, 23.] , Ederth 0, Chan et al. 0, and 



Bressi et al. |26| (cf. also the recent review paper of Lambrecht and Reynaud [27|), represents 
important progress in this field. Especially the experiment of Bressi et al. is of interest in the 
present context, since it deals directly with the Casimir force between metal surfaces that are 
parallel, and so avoids use of complicating factors such as the proximity force theorem p8l |- 
which nevertheless seems well understood. This experiment is fraught with experimental 
difficulties (related to keeping the plates sufficiently parallel), so the accuracy is claimed by 
the authors to be moderate (15%), but it is to be hoped that this accuracy will soon be 
improved. Several other related papers have appeared recently, discussing the interpretation 
of the mentioned experiments as well as more general aspects of finite temperature Casimir 
theory H H H HI H H Q . 

Our purpose in the present paper is to analyze the Casimir temperature problem anew, 
assuming conventional parallel-plate geometry from the outset, therewith avoiding the spher- 
ical Bessel functions that become necessary if spherical geometry is contemplated. In par- 
ticular, we will focus attention on the m = TE mode. Let us summarize our results: 

It is useful to distinguish between two different classes of metals. The first class, which 
we will call "ideal" metals, is characterized by a refractive index n = y/e = oo for all 
frequencies. It implies that the refiection coefficient r2 mentioned above is unity for all 

This corresponds to the traditional recipe ^ and |H above when handling the two-limit 
problem for metals. It means that the m = TE mode contributes to the Casimir force just 
the same amount as does the transverse magnetic (TM) mode. 

The obvious drawback of this "ideal" metal is that it does not occur in nature. And this 
brings us to the second class, which is the one of real metals, in which case we must observe an 
appropriate dispersion relation, especially at low frequencies. It is most commonly assumed 
that the most appropriate dispersion when — * is the Drude relation, Eq. (jl.lj) . As we 
will show, the Drude model implies that the m = TE mode does not contribute. The 
total m = free energy for a real metal becomes accordingly one half of the conventional 
expression. In contradistinction to recent statements in the literature [1, 0, 13 we find that 
there exists no physical difficulty or ambiguity associated with the vanishing coefficient r2 
at ^ = 0. This is so because r2 goes to zero smoothly when ^ — 0, as long as the transverse 
wave vector is nonvanishing. (If k_|_ is precisely zero, there occurs a singularity in the 
reflection coefficient, but this has no physical importance since this point is of measure zero 
in the integral over k^.) Our present results are in agreement with Refs. Esj]. as well as 
with Bostrom and Sernelius |l3|. 

A different view has recently been put forward by Torgerson and Lamoreaux [s^- They 
argue that the Drude-model behavior does not accurately represent the TE zero mode, 
which necessarily has a vanishing tangential component at the surface of a perfect conduc- 
tor. They point to the necessity of taking the flnite thickness of the metallic coatings into 
account. Their arguments seem to imply that the conventional temperature dependence is 
correct. However, in our opinion electrostatic considerations of this kind do not solve the 
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zero temperature problem; what is required to incorporate temperature dependence is an 
analytic continuation into imaginary frequencies of Green's functions referring to nonzero 
wavenumber. 

Before embarking on the calculations let us emphasize the following point: The occurrence 
of the m = mode only once instead of twice is understandable physically. This mode is 
precisely the TM static mode, corresponding to the electric field being perpendicular to the 
two metal plates. It is the natural ground-state mode present when ( = 0. Actually, in 



Sec. Ill of Ref. [12| we showed how the uniqueness of the static mode emerges naturally, 
using statistical mechanical considerations. 

The outline of our paper is the following. In the next section we show why the exclusion 
of the TE zero mode seems to lead to an observable temperature correction to the force 
between real metal plates, and worse, seems to imply a violation of the third law of thermo- 
dynamics. In Sec. mil we expand on the situation of an "ideal" metal in the sense described 
above, and calculate the Casimir free energy, entropy, and internal energy via a somewhat 
unconventional route. Equivalence with earlier results is demonstrated. In Sec. IIVI we in- 
troduce a new and simplified model to illustrate the Casimir problem, based essentially on 
statistical mechanics. In this model the system is replaced by two harmonic oscillators (the 
two media) that interact via a third oscillator (the electromagnetic field). Depending upon 
the form of the interaction we then have two situations. The first is the one where the 
induced interaction (or free energy), which is negative, increases linearly in magnitude with 
temperature in the classical limit. The other situation, which is more unexpected, is where 
the induced interaction vanishes in the classical limit. These two situations can be regarded 
as analogous to the behavior of the TM and TE modes. We also consider a strongly sim- 
plified case of real metals, and show how in such a case the contribution to entropy goes to 
zero smoothly as T — >■ 0. Arguing on basis of the Euler-Maclaurin formula we find this to be 
a general property (except in the idealized metal limit). We then go on to present numerical 
results based on the dispersion relation for gold, and obtain results qualitatively in accord 
with our analytical model. In the Appendices the smoothness of the reflection coefficient 
r2, and of the TE Green's function, in the limit C ^ is explicitly demonstrated. We also 
discuss the temperature dependence of the relaxation frequency, z^(T). We conclude that a 
linear temperature dependence should be observable in room temperature experiments. 

In this paper we use natural units, h = c = ks = 1- 



II. TEMPERATURE EFFECT FOR METAL PLATES 

We begin by reviewing how temperature effects are incorporated into the expression 
for the force between parallel dielectric (or conducting) plates separated by a distance a. 
To obtain the finite temperature Casimir force from the zero-temperature expression, one 
conventionally makes the following substitution in the imaginary frequency, 

27rm 

C^Cm = (2.1a) 
and replaces the integral over frequencies by a sum, 
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This reflects the requirement that thermal Green's functions be periodic in imaginary time 
with period P fsT^]. Suppose we write the finite-temperature force/area as [for the exphcit 
form, see Eq. ()3.1|) below] 



oo 



= E (2.2) 

where the prime on the summation sign means that the m = term is counted with half 
weight. To get the low temperature limit, one can use the Euler-Maclaurin (EM) sum 
formula, 

J2fik)= f{k)dk + -f{0)-J277^/''-'\0), (2.3) 

k=0 ^ q=l ^^^>- 

where Bn is the nth Bernoulli number. This means here, with half-weight for the m = 
term, 

/■oo -1 -1 °° r) 

= I f{m) dm - ^/(O) + ^/(O) - J2 (iy!/^''"'^(0)- (^.4) 

It is noteworthy that the terms involving /(O) cancel in Eq. ()2.4|1 . The reason for this is 
that the EM formula equates an integral to its trapezoidal-rule approximation plus a series 
of corrections; thus the 1/2 for m = in Eq. ()2.2p is built in automatically. For a perfect 
conductor 

Of course, the integral in Eq. ()2.4j) is just the inverse of the finite-temperature prescription 
dnO), and gives the zero-temperature result. The only nonzero odd derivative occurring is 

r(0) = (2.6) 

which gives a Stefan's law type of term, seen in Eq. (j2.1(J|) below. 

The problem is that the EM formula only applies if f{m) is continuous. If we follow 
the argument of Ref. [3, 3, and take the ei,2 — * C)0 limit at the end (ei,2 are the 



permittivities of the two parallel dielectric slabs), this is not the case, and for the TE mode 

/o = 0, (2.7a) 

Then we have to modify the argument as follows: 

oo oo 



m=0 


— ^ frn 

m=l 


oo 
m=0 





(2.8) 



where fm is defined by continuity. 



fm=< r ^"^ . ^ n (2.9) 



Then by using the EM formula, 



27r 7o 87r/3a3 45 \I3 

C(3) 



240a4 



16 / ^ 



T, aT < 1. (2.10) 



The same result for the low-temperature limit is extracted through use of the Poisson sum 
formula, as, for example, discussed in Ref. |l|. Let us refer to these results, with the TE 
zero mode excluded, as the modified ideal metal model. 

Exclusion of the TE zero mode will reduce the linear dependence at high temperature 
by a factor of two, but this is not observable by present experiments. The main problem, 
however, is that it adds a linear term at low temperature, which is given in Eq. ()2.10p . up 
to exponentially small corrections 0. 

There are apparently two serious problems with the result (j2.1(jp : 

• It would seem to be ruled out by experiment. The ratio of the linear term to the T = 
term is 

A = ^^^^aT = l.maT, (2.11a) 
or putting in the numbers (300 K = (38.7)^^ eV, he = 197 MeV fm) 

A = 0-15f;3^VAV (2-llb) 



300 K J \lfim 

or as Klimchitskaya observed j3^, there is a 15% effect at room temperature at a 
separation of one micron. One would have expected this to have been been seen by 
Lamoreaux ; his experiment was reported to be in agreement with the conventional 
theoretical prediction at the level of 5%. 

Another serious problem is the apparent thermodynamic inconsistency. A linear term 
in the force implies a linear term in the free energy (per unit area). 



^o + S^T, aT«l, (2.12) 



which implies a nonzero contribution to the entropy/area at zero temperature: 

Taken at face value, this statement appears to be incorrect. We will discuss this problem 
more closely in Sec. I1V[ and will find that although a linear temperature dependence will 
occur at room temperature, the entropy will go to zero as the temperature goes to zero. 
The point is that the free energy F for a finite e always will have a zero slope at T = 0, 
thus ensuring that S' = at T = 0. The apparent conflict with Eq. ()2.13j) or Eq. ()2.1()|1 
is due to the fact that the curvature of F[T) near T = becomes infinite when e — > oo. 
So Eqs. ()2.12p and ()2.13p . corresponding to the modified ideal metal model, describe real 
metals approximately only for low, but not zero temperature — See, for example, Eq. ()4.14p . 
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III. CASIMIR FREE ENERGY, ENTROPY, AND INTERNAL ENERGY 



The Casimir surface force density between two dielectric plates separated by a distance 
a can be written as 



T °° poo 



4 „-2qa f) p-2qa 



1 - AmC-^i" 1 - Bme-^i" 



(3.1) 



(We follow the conventions of Ref. [39| and further references therein; here we further set 
h = c = 1.) The relation between q and the transverse wave vector kj^ is = k'j_ + 
where (m = 27rm//3. Furthermore 

A„=(fP^)\ B„=f^)'. (3.2a) 

\ep + sj \s+pj 

= e-l+p\ P=j-, (3.2b) 

with e{i(m) being the permittivity. Note that whenever e is constant, the and Bm depend 
on m and q only in the combination p, 

AUq)=A{p), B^{q) = B{p). (3.3) 



(This result may also be found in standard references such as Ref. [l|.) 

The free energy F per unit area can be obtained from Eq. ()3.H) by integration with 
respect to a since JF^ = —dF/da. We get ll^ 



1 _ _ 

PF = -J2' Ml-X™) + ln{l-X^^)]qdq, (3.4a) 

where 

A™ = A^e-2'?^ A^^^ = Bme-''i\ (3.4b) 

(In the notation of Ref. [3, A, = A™, A = A™.) 

From thermodynamics the entropy 5* and internal energy U (both per unit area) are 
related to F hj F = U — TS, implying 

S = -^, and thus U=^^. (3.5) 

As mentioned above the behaviour of as T — has been disputed, especially for metals 
where e oo. We now see the mathematical root of the problem: The quantities Am = 
Bm — 1 in the e — oo limit except that -Bq = for any finite e. So the question has been 
whether i?o = or i?o = 1 or something in between should be used in this limit as results 
will differ for finite T, producing, as we saw above, a difference in the force linear in T. The 
corresponding difference in entropy will thus be nonzero. Such a difference would lead to 
a violation of the third law of thermodynamics, which states that the entropy of a system 
with a nondegenerate ground state should be zero at T = 0. Inclusion of the interaction 
between the plates at different separations cannot change this general property. We will 
show that this discrepancy vanishes when the limit e — oo is considered carefully, by using 
the Euler-Maclaurin summation formula. Also, we will perform explicit analytic evaluation 
for any T for metallic plates in the case where e ^ oo for all C,. 

We will consider this latter case first. It is the case of "ideal" metals mentioned in Sec. H] 
and already considered briefly in Sec. |Tl| 
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A. "Ideal" metals 



With e = oo we have Am = Bm = 1 where we now also put Bq = 1, i.e., A™ = A"""^ = 
g-2qa rji^ remove the ^-dependence in the lower limit of integration in Eq. ()3.4ap . it is 
convenient to use the quantity p of Eq. ()3.2b|) as a new variable. Expanding the logarithmic 
terms in Eq. (j3.4aj) and keeping only the leading term, we get the task of calculating 



1 °° poo 



dp, (3.6) 



m=0 

where 



Carrying out the integration in Eq. (j3.f)j) we obtain 



27ra 



^ oo 

^1 = 7^2^'^^' (3.8a) 



(2a) 
with 

Lm = (27m + l)e-2T"^. (3.8b) 

(It is easy to check that this result is correct at m = 0, where p is not defined.) We encounter 
the following sums 

00 

so(^) = 2^'e-2^™ = coth7, (3.9a) 



m=0 
00 



00 r\ji 

2j](27m)V^-^ = (-7)^|;^, (3.9b) 



m=0 



SO that 



7 

si = ——n—, (3.10a) 
smh 7 

272 cosh 7 /QinM 

S2 = — —3 , (3.10b) 

smh 7 



36 + 4 sinh^ 7 1 n ^ 

S3 = 7 -Ti • (3.10c) 

smh 7 

The quantity Ii is given by the first two of these sums, 

hil) = j^[siil) + soii)]. (3.11) 

Alternatively, one could just first perform the summation in Eq. ()3.6p (for m > 1) and 
then integrate. This summation yields 52(7^). By subsequently integrating S2 by parts the 
quantity (si + sq) in Eq. (j3.11|) is recovered (adding the m = term separately). 
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By further expansion of the logarithm in Eq. ()3.4ap one obtains terms \^/k to be inte- 
grated and summed hke Eq. ()3.6p . Performing the same steps as before, we find that the 
result ()3.11|) generalizes to 

1 °° 1 

^=-8^5:pN.W + ^«(7*)]. (3.12) 

^ k=l 

valid for arbitrary temperature. 

The surface force per area p.ip can now be obtained via = —dF/ da utilizing 7 oc a 
[Eq. dnHI)]. This yields 

00 

= E ^[^2(7^) + 2s,{^k) + 2.o(7fc)]. (3.13) 

k=l 

The same result is also obtained by evaluating expression (|3.1|) (with = -Bm = 1) in 
the same way as expression (j3.4ap for F was evaluated above. Using the second method, 
mentioned below Eq. ()3.11|) . one finds that the integration of 53(7^) /p yields the combination 
of Si present in Eq. ()3.13p . 

Considering the T limit, which implies the 7 — > limit, one obtains 

= = (314) 

^ k=l ' 

using the limiting values of expressions (j3.9aj) . (j3.1(Ja|) . and ()3.10b|) . This is the well known 
Casimir result for idealized metallic plates at T = 0, seen in Eq. ()2.10j) . 

The internal energy U is now found from Eqs. ()3.5p . ()3.7p . and ()3.12|) to be 

and similarly an expression for the entropy S follows from 

TJ p 1 00^ 

S = -2na— = ^ = E p[^2 W - ^1 W - -^oilk)]. (3.16) 

' k=l 

with Eqs. flTT^ and (ITTHIl inserted. 

Now we can analyze the thermodynamic quantities in the low temperature limit using 
the properties of as defined by Eqs. (j3.9aj) - (j3.1(Jc|) . We have for low temperature,^ where 
7 oc T ^ 

1 1 





1 




So = 




+ 




7 






1 




Si = 






7 






2 




S2 = 








7 






6 




S3 = 




+ 




7 





_ 0,3 

3 ' 45 ' 

1 1 3 

-7 H 7 

3 ' 15 ' 

2 

^y-.... (3.17d) 



(3.17a) 
(3.17b) 

7^ + ..., (3.17c) 



^ Actually, for a room-temperature experiment, 7 need not be small. For T — 300 K and a — l^m., 
7 = 0.823. 
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Inserting this into expressions ()3.12|) . ()3.13|) or ()3.15|) one finds that the terms hnear in 7 
vanish.^ Thus the entropy ()3.16p vanishes, as it should in accordance with the third law of 
thermodynamics. 

To obtain the leading correction to the T = result for finite T one must consider the 7^ 
term in the power series expansion of the summand in Eq. (|3.13p . However, the summation 
of this term with respect to k diverges,^ because the expansion of Snipfk) is not valid for 
large k. For small 7 one can instead integrate, without expanding, using the Euler-Maclaurin 
summation formula ()2.3|) to obtain a finite correction to the zero-temperature result. Using 
Eq. ()2.3|) to evaluate expression ()3.13p . the 7^0 expression ()3.14p has to be subtracted to 
make /(O) finite. Putting x = jk we have, apart from a prefactor. 



X-' 



6 

S2(x) + 2si(x) + 2so(x) 

X 



(3.18) 



with /(O) = —2/45 in view of the expansions ()3.17aj) - ()3.17c|l . Integrating and using expres- 
sions (jTT?^ . (ITTnill . (I3l0b|) . we obtain 



f{x) dx = — - 



'0 ^ 
Including the T = result ()3.14|1 we thus find 

1 



siix) + soix) - 



X 



(3.19) 



240a4 



7 90 2 



^"5 7 



aT < 1, 



(3.20) 



where we have inserted expression ()3.7|) for 7 and noted that there is no /c = term in 
Eq. ()3.13|) . i.e., /(O) is to be subtracted from expression ()2.3|1 . All the odd derivatives 
in the Euler-Maclaurin formula vanish because f{x) is even. It should be noted that the 
expression for J-'^ is in agreement with what has been found earlier [cf. Eq. (I2.1()|l ]. via 
alternative methods, by Milton jl|, Klimchitskaya and Mostepanenko 0], Sauer |40j, Mehra 
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and others, where the exponentially small correction to the above formula is also given. 
The free energy ()3.12|) can be obtained from JF^ = —dF/da, but this leaves a temperature 
dependent constant of integration. So instead we make use of the method above, where from 
Eq. (HTT^ 



x^ 



sAx) + Sq(x) 

X 



(3.21) 



and where now /(O) = 2/45. With Eq. 1)3.211) we get a nonzero integral 



C 



f{x) dx 



x-^ 



X 

3 



°° 1 d 

'tic doc 

cothx I dx 



1 , 1 

— coth X 

X 3 



x^ 



dx 



(3.22) 



^ This is actually stronger than necessary to insure vanishing entropy, since such terms would give terms 
in the energy or free energy. 

For this reason, the alternate expression (3.35) in Ref. |l| might be preferred. See Eq. H3.39|l below. 



10 



using partial integration. The integral ()3.22|) may be easily evaluated by contour methods. 
Due to symmetry the integral can be extended to minus infinity and then the contour of 
integration can be distorted into one which encircles the poles along the positive imaginary 
axis. Since coth z has poles a.t z = mm with m integer we get^ 

1 °° -1 1 
C = -2m V —— = -C(3). (3.23) 

m=l ^ ^ 

In view of this result as well as Eq. ()3.14|) we obtain for the free energy [dk = dx/j) 

1 /2 7r4 3 1 



7r/5a2 V79O ' V7 2 



+ 7M--o/(o) 



l + 45('^V^-f^Y l> -T«l. (3.24) 



720a3 \' \(3 J \ (3 J 

This result, including its exponentially small correction, is given in Ref. |lj and references 
therein. The internal energy U, which can be most easily be evaluated using Eq. (j3.5p . can 
also be computed by the method above, starting from the sum (j3.15|) . Then 

/(-) = i (».(-) - 1) = 4^ + - • f^-^^) 

with /(O) = —2/15. Partial integration replaces the C of Eq. ()3.22|) with —2C, and we 
obtain 



U 



7r2 



aT < 1. (3.26) 



720a3 

With Eq. ()3.16|) the entropy thus becomes (recall that i?o = 1 is assumed) 

U-F 3C(3) 2 47r2a 

B. Equivalence with earlier results 

Equivalence with previous derivations can be shown for any 7. It is then convenient to 
utilize the Poisson summation formula. If c{k) is the Fourier transform of c(x), defined by 

/oo 
dxc{x)e'''\ (3.28) 
-00 

then 

00 00 

c{n) = d{2nm). (3.29) 

?i=— 00 m=— 00 



^ This low temperature dependence in F, which does not contribute to the force, is determined by the 
linear high temperature behavior of — see Ref. Q^, Sec. 3.2.1. 
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With c{x) = e~^'''l^l one finds 

c(27rm) 

Thus 



^-2j\x\+2iTmxi 



7 



7^ + (vrm)^ 



m=—oo 



7 



7^ + (vrm)" 



^ e-2^l"l = coth7, 



(3.30) 



(3.31) 



the famihar cotangent expansion, which can be verified in many different ways (cf. Ref. |42l|). 
In Eqs. fl3.12|) and ()3.13|) one of the sums is [so{x) = cothx] 



oo oo 



•S'o = ^ p coth{'yk) = Somk, 



k=l 



k=l m=—oo 



where with Eq. ()3.31|1 



5-, 



7/c 



Omk 



(3.32a) 



u = 7r^m/7. (3.32b) 



Summation first with respect to k where also the result ()3.3H) is utihzed then gives 



5, 



Om 



fc=i 



Omk 



mu 



71^ vrV , 1 

cotli u 

6 2m \ u 



(3.33) 



In the limit 7^0 only the m = term remains, and we get the T = result if we use the 
expansion ()3.17a|l {u 0) 



5', 



00 



mu 



2^7 V 457 907' 



(3.34) 



which is consistent with the l/fc^ sum occurring in Eq. ()3.14|1 . 

To obtain the free energy F and the force JF^ there are sums 5*1 and S2 that follow from 
the Si and S2 of Eqs. (|3.1(Ja|) and ()3.10b|) . And like Eqs. (j3.32aj) the relations between the 
various Si lead to 

S ^ 



where g{u) = Som- Also: 



-7 -^9 = ug , 
07 



So to obtain JF^ we need, because 



d_ 
du 



2 + 4m— + u'— 1 —g{u) = g"{u), 



du^ ) u 



(3.35) 



(3.36) 



(3.37) 



the combination 

S2m + 25*1^ + 2S'o„ 



6m du^ 
/ 1 



u — 3(cothM 



u 

cosh-u \ m^o vr^ M vr"^ 1 



m \u-^ sinh M 



m 15 15 7 



(3.38) 
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Altogether, restricting m to positive values due to symmetry, the expression ()3.13p can be 
reexpressed as {u = n^m/'y, 7 = 2na/P) 



IT 



m=l ^ 



coshw 
u sinh^ u 



(3.39) 



which is the desired known expression. (For example, compare Eq. (3.35) of Ref. [l|.) 
To calculate the free energy ()3.12|) one likewise needs 



Sim + Si, 



Om 



6m du 



1-3 



71 

2m 
2 



coth u 



+ 



f cothw 


^2) 


V u 




1 


2 


u sinh^ u 





2m V 45 15 



u 



!L1 
457' 



(3.40) 



Thus the free energy becomes 

^2 



7r 



720a3 



1 + 45 5:1 



/ coth 



u 



+ 



m=l 



V? sinh u 



(3.41) 



Compared with the small T or 7 expansion ()3.24j) it is clear that the last term of Eq. ()3.41|) 
gives the = j3~^ term of ()3.24|) . The coefficient C can also be identified from Eq. ()3.4H) . 
As coth-u 1 when 7 — we must have, when comparing with Eq. ()3.24|) . 



7\4 45C 



TT 



7 



00 ^ 



m=l 



or 



C 



45 



rC(3), 



7 



, 00 -I 

Z / jyj-^ 

m=l 



(3.42a) 



(3.42b) 



which is in agreement with Eq. (j3.23j) . 



IV. FINITE PERMITTIVITY. REAL METALS 
A. Two harmonic oscillator models 

With finite permittivity e the and of Eq. ()3.2a|) will vary with p. Especially 
Bm — > as p 00 or — > (Cm = 27rm//9). In the high temperature or classical limit 
only the Matsubara frequency C = (or m = 0) can contribute as /3 — >■ 0. Thus, in the 
classical limit one has the result that the TE mode does not contribute at all. Physically, this 
means that the temperature becomes so high that only the static dipole-dipole interaction 
contributes (the C ^ limit of the TM mode). In our opinion this somewhat unexpected 
behaviour is related to the peculiar type of interaction that exists between the canonical 
momentum p of a particle and the electromagnetic vector potential A(r, t), which for a 
particle of mass m and charge g is (p — gA)2/2m. In addition to the standard cross term 
interaction p ■ A this also implies an interaction A^. 
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As an illustration of the above we can consider two models, in each of which two harmonic 
oscillators interact via a third one. These oscillators represent a simplified picture of our 
polarizable parallel plates interacting via the electromagnetic field. The classical partition 
function of a harmonic oscillator with frequency uj is const/ (/Jcj) ~ l/y/u^, which gives a free 
energy ~ hi{uj'^). Thus for three noninteracting harmonic oscillators the inverse partition 
function is proportional to y^, where 

Q = aia2a3, (4.1a) 

with 

a, = u;l (2 = 1,2,3). (4.1b) 

(The quantity corresponds to kj_ above.) By quantization using the path integral method 
(4^ 43l |. the classical system is split into a set of harmonic oscillator systems described by 



Matsubara frequencies. Expression ()4.1ap is replaced by 

Q = AiA^As, (4.2a) 

where 

(For real frequencies, uj = i(, 1/Ai determines the response to an external oscillating force 
acting on the oscillator.) 

Now add interactions, of strength proportional to c, between the third oscillator and the 
other two. The usual form of this interaction is cXiXj, where Xi and Xj are coordinates. Let 
this constitute the first model, which is analogous to the TM mode. Then the quantity Q 
becomes the determinant of the matrix, 



A1A2A3 - c\Ai + A 



where 



Ai c 
Q = A2 c 
c c As 

A,A,Ml-DMl-D,)(l-^^-^g-^), (4.3a) 

A = (. = 1.2). (4.3b) 

The quantum free energy for this system of three coupled oscillators is given by summing 
over the Matsubara frequencies, as in Eq. ()3.4a|) : 

1 ^ 

f3F=-\im J2(}nQ{U) + 3W), (4.4) 

m=l 

where rj = j3/N and ("^ is replaced by 2(1 — cos{(r]))/r]'^ (= (^^ + ...) in the A1A2A3 term 
of Eq. ()4.3a|l . The limiting procedure ^ 00 is required to make the full free energy 
well defined. This means that the path integral representation of a harmonic oscillator is 
discretized by dividing the imaginary time of periodicity (3 into pieces each of length rj as 
done in Ref. |42|. There, in an appendix an explicit evaluation was performed for one single 
oscillator. 
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The various factors in Eq. ()4.3a|) can be interpreted as follows: The product A1A2A3 
corresponds to the noninteracting system, the next two factors represent the result of in- 
teraction of single oscillators with the third one, while the last one is the contribution from 
the induced interaction between the two single oscillators via the third one. The logarithm 
of the last term is the analogue of the Casimir free energy. In this respect the term c^/A^ 
represents the induced interaction. Furthermore the 1/aj {i = 1,2) represents the "bare" 
polarizability of noninteracting particles which for nonzero ( becomes 1/Ai. Due to interac- 
tion with the "radiation" field this polarizability is modified into l/{Ai{l — Di)) {i = 1,2), 
where Di represents a "radiation" reaction from the "field" upon each single oscillator. 

The above represents the ordinary situation, analogous to the TM mode. To model the 
TE mode, we can consider an analogy with the electromagnetic interaction in which the 
third oscillator can interact with the momenta of the first two. The analogous interaction 
will be {pi — const. Xs)^ /2m j {i = 1,2; rrij is mass), including the unperturbed p] term. 
By evaluation of the classical partition function one now finds that the interaction from 
const. X3 has no infiuence. (This is the analogue of classical diamagnetism which is equal to 
zero, as const. 0:3 is seen to have no infiuence on the result when pi is integrated first.) 

Quantum mechanically, the problem is a bit more complex. However, we can now ex- 
change the roles of momenta and coordinates of the first two oscillators, i.e., we introduce a 
momentum representation. Then the interaction with the third oscillator can be written as 
(z = l,2) 

\ 2 / 2 



const, tti { Xi 0:3 = const. I a^Xj — 2cXiXs H 0:3 . (4.5) 

Now the last quadratic term adds to the energy of the third oscillator alone. Thus, compared 
to the first model considered above, as is changed while the other at remain unchanged: 

^ ^3 + c^/fli + c^/a2. (4.6a) 

Likewise in the quantum case 

A3 A3 + cVai + c^/a2. (4.6b) 

The quantity Q can still be written in the form ()4.3a|) . but due to the change of 03, the 
(l/Ai) {i = 1,2) is replaced by l/Ai — l/a, = — C^/(ajAj) when evaluating D^, i.e., 

= 'a.ia. + OA', ^^•'^^ 

The induced (analogous to the Casimir) free energy is again given by the logarithm of the 
third term in Eq. ()4.3a|) . At zero and finite temperatures the latter logarithm is negative, 
and the free energy 

T ^ , f D1D2 \ , , 

- V n 1 - 4.8 

is negative. Note that here the limiting procedure of Eq. ()4.4j) is not needed as sums for 
free energy differences converge, without difficulties. In the classical limit, however, the 
induced free energy becomes equal to zero [Di implies that we get the logarithm of 
unity). We note the analogy: At high temperatures the same is true for the TE mode in the 
Casimir effect. There exists thus at least somewhere a finite temperature interval for which 
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the Casimir free energy increases with increasing temperature. In turn, this means that the 
Casimir entropy S = —OF/dT becomes negative in this interval. 

This is a counterintuitive effect, but is physically due to the fact that we are dealing with 
the induced interaction part of the free energy of a composite system. We cannot apply usual 
thermodynamic restrictions such as positiveness of entropy to a "subsystem" of this sort. 
There exists actually a striking analogy with the peculiar formal properties one encounters 
in connection with the theory of the electromagnetic field in a continuous medium. The 
electromagnetic energy-momentum tensor that experimentally turns out to be definitely 
the best alternative when dealing with high-frequency effects, is the Minkowski tensor (cf., 
for instance, Ref. This tensor is however nonsymmetric, apparently breaking general 

conservation principles for angular momentum. The reason why this peculiar behaviour is 
yet quite legitimate physically, is that phenomenological electro dynamic theory is dealing 
only with a subsystem (the field itself plus its interaction with matter), and we cannot apply 
the same formal restrictions on it as we could if the system were closed. 



B. Real metal 

In the limit of an ideal metal {e —>■ oo) the traditional (SDM) prescription, as mentioned 
in the Introduction, implies that Am = = 1 for all m. In addition, as also mentioned 
previously, thermodynamic arguments have been given, claiming that the entropy does not 
become zero at T = in violation of the third law of thermodynamics if i?o = is used js^ . 
However, we do not find this to be the case; as we will show below, the entropy will be zero 
as required at T = 0, even for a metal that is not idealized and where one bases the analysis 
on the value So = 0. 

Let us go back to Eq. ()3.6|) . That equation was obtained by expanding Eq. ()3.4a|) to first 
order in A under the assumption that Am = Bm = 1- Doing the same expansion for finite 
permittivity, we obtain an integrand which contains a term with a factor B^, (or A^) that 
varies with p = q/Cm such that B^ when p ^ oo. Expanding Eq. (j3.4aj) to higher 
order one obtains likewise powers of Bm which, because Bm < 1, become less important as 
compared to the case of an ideal metal (where Bm = 1). One can first consider the case 
where e is independent of (. When e is large one can use as a rough approximation 

fl, p<V£, 

" \ 0, p > v^. ^^-^^ 

This simple expression for Bm is intended to show essential features that will be obtained 
more accurately in a detailed numerical calculation. With this, Eq. ()3.9ap (neglecting the 
influence of Am) will turn into 

Soil) Soil) - soiy^-f) = coth7 - coth7c, (4.10a) 

with similar modifications for Si (z = 1,2,3). Here 

Ic = 7v^ (4.10b) 

is an effective sharp cutoff limit for the integral, a crude model for what should be a gradual 
cutoff for the integral of interest. [A gradual cutoff will only modify the last term of ()4.10a|l 
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into a sum or integral over terms with varying 7c. Namely, with varying B = B{p), Eq. ()3.6p . 
if we recall the comment below Eq. ()3.1H) . changes into (5(1) ~ 1 for e large) 

^1 = 7W~^ [ B{p)s2{np) — 
{2ay Ji p 

1 1 /"°° 

= (2^ + + (2^ / + si(7p)]S'(p) dp, (4.11) 

using partial integration. The approximation (j4.9p means that B'{p) = —6{p — y/e).] 

As we did to obtain Eq. ()3.27|1 . we carry out the sum over k in Eq. (jH.lfij) while assuming 
e sufficiently large such that approximation ()4.9|1 can be used. Then as in Eq. ()4.1()aj) one 
obtains the previous result minus a term with 7 — > 7c. Keeping only the leading term, 
Eq. ()3.27|) is modified into 

[However, to be more accurate Bm = ((v^ — l)/(v^ + 1))^ for p = 1 and thus Bm < 1 for 
p < y/e. When this is taken into account, we find that 5*™ oc — ae^/^T^ in a more narrow 
region, e^^^ aT < 1, but that Eq. KT^ holds for e'^^^ <^ aT <^ e^^/^ ] 

Thus the entropy approaches zero as the temperature goes to zero. As e increases the T- 
dependence becomes more singular, because the region in which Eq. ()4.12p is valid becomes 
more and more narrow, but the value at T = stays fixed at zero also in the limit e — > cxo. 
This contrasts the ideal metal result ()2.13|) where e = cxd is used. 

Again, we note the counterintuitive negative contribution from the TE mode. As men- 
tioned earlier, this does not violate the laws of thermodynamics and can be understood in 
terms of the oscillator model analysed in some detail in Sec. II V Al Only the total entropy 
has to increase with increasing temperature. And this is the case for the inverse partition 
function ()4.3ap which represents three interacting harmonic oscillators where the Di are 
given by Eq. ()4.7|) . Although the induced entropy becomes negative at least in some finite 
temperature region the total entropy will behave properly, as the total system can be de- 
composed into three independent harmonic oscillators represented by the eigenvalues of the 
matrix ()4.3ap with Ai replaced by Oj {i = 1,2), and furthermore A^, replaced by the right 
hand side of Eq. (|4.6ap . 

With the simplification (|4.9p for the TE-mode the free energy can be easily expressed in 
terms of the "ideal" metal case analysed in Sec. 1111 Al Let the "ideal" metal free energy be 
F = Fi{T). From Eq. ()3.7|) 7 oc T. Now the magnification of 7 to 7c as in Eq. ()4.10b|) and 
insertion of it in Eq. ()3.12|) will change the corresponding free energy to {'j /'-fc)Fj{T'-fc/'j) = 
FI[^/eT)/ ^/6. The TM- and TE-modes both contribute the same amounts to ()3.12p . Thus 
with Eq. ()4.9|) the free energy will be 

F = F{T) = Fj{T) - ^Fji^T). (4.13) 

From this we have (keeping in each case only the leading temperature correction) 

1 - l/(2v^))F,(0) - ^(2 - e)T\ < aT « 1/^, 
F(T) = { Fi{0) + KiT/2, 1/v^ < aT < 1, (4.14) 

-KjT/2, 1 < aT. 
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where the constant Ki = ({3) / (Sna^) is the magnitude of the slope of the hnear dependence 
of the high temperature resuh of the the "ideal" metal {Fi{0) = —7^5^ < 0). Thus for 
high temperatures non-ideal or realistic metals yield one half of the "ideal" metal result. 
The intermediate form, which holds at room temperature, is the same as seen in Eq. ()2.12p . 
Again, we see that in the y/eaT <^ 1 regime the result ()4.12|) for the entropy holds. [Equation 
(|4.14p includes the TM mode as well.] 

Now, e usually depends on (. But this will not change our conclusions from Eq. (|4.12p . To 
see this we can go back to expression which followed from expansion of the logarithmic 
term in the free energy ()3.4a|) . In the general case, the coefficients Am and Bm, which are less 
than 1, should be included in Eq. ()3.6p . and powers of them will occur in the evaluation of the 
terms contributing to the free energy for k > 1. These factors will all be smooth functions of 
( except for the case of an idealized metal where becomes discontinuous aX ( = 0. This 
smoothness is also valid for the Drude formula discussed in Appendix^ With and B^, 
included, Eq. (I3.4a|) can be summed with respect to Cm = ^nm/P, and the Euler-Maclaurin 
formula (j2.3|) can again be applied. (Equation (j3.6p with B^ included is not applicable in 
this situation as we remarked there because C is of relevance here.) If e stays finite 
when ( —>■ the result clearly will be the same as that given above. However, for a real 
metal where e — 00 as C the situation is more subtle. For the case of an ideal metal 
considered in Sec. |nj the first derivative /'(O) was zero while /'"(O) of Eq. ()2.6|) was nonzero. 
By similar application of the Euler-Maclaurin formula to the free energy ()3.4a|) instead of 
the force p.ip . the same will be true. For a real metal obeying the Drude dispersion relation 
(with u 0) the first derivative /'(O) continues to be zero due to the ( dependence 
of Bm, Bm ~ Cm; according to Eq. ()A4|) . Thus, quite generally, we expect a (or T^) 
correction to the free energy at sufficiently low temperature. 



C. Gold as a numerical example 

Let us go back to Eq. (j3.H) for the surface force density, making use of the best available 
experimental results for e{i() as input when calculating the coefficients Am and Bm- We 
choose gold as an example. Useful information about the real and imaginary parts, n' and 
n", of the complex permittivity n = n' + in", versus the real frequency u, is given in Palik's 
book [4^ and similar sources. The range of photon energies given in Ref. |43] is from 0.1 
eV to 10^ eV. (The conversion factor 

1 eV = 1.519 X 10^^ rad/s (4.15) 

is useful to have in mind.) When n' and n" are known the permittivity e{iQ along the 
positive imaginary frequency axis, which is a real quantity, can be calculated by means of 
the Kramers- Kronig relations. 

Figure m shows how 6{iQ varies with ( over seven decades, ( G [10^^,10^®] rad/s. The 
curve was given in an earlier paper ji^], and is reproduced here for convenience. (We are 
grateful to A. Lambrecht and S. Reynaud for having given us the results of their accurate 
calculations.) At low photon energies, below about 1 eV, the data are well described by the 
Drude model, Eq. ()l.ip . in which the input parameters have the values [s^l 

Up = 9.0eV, z/ = 35 meV. (4.16) 

These values refer to room temperature. The curve in Fig. ^ shows a monotonic decrease of 
e{iQ with increasing (, as any permittivity along the positive imaginary axis has to follow 
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10" 10" 10" 10" 10" 10'° 10" 10'° 10" 

C (rad/s) 

FIG. 1: Full line: Permittivity e(iC) ^ function of imaginary frequency (" for gold. The curve is 
calculated on the basis of experimental data. Courtesy of Astrid Lambrecht and Serge Reynaud. 
Broken lines: e(«C) versus (" with T as parameter, based upon the temperature dependent Drude 
model; cf. Appendix^ The upper curve is for T = 10 K; the lower is for T = 300 K, which for 
energies below 1 eV (1.5 x 10^^ rad/s) nicely fits the experimental data. Both curves are below the 
experimental one for C > 2 X 10^^ rad/s. 

according to thermodynamical requirements. The two broken curves in the figure show, 
for comparison, how 6{i(,T) varies with frequency if we accept the Drude model for all 
frequencies, and include the temperature dependence of the relaxation frequency with T as 
a parameter. Cf. Appendix O For T = 300 K, the Drude curve is seen to be good for all 
frequencies up to C ~ 2 x 10^^ rad/s; for higher ( it gives too low values of e. Both Drude 
curves, for T = 10 K and T = 300 K, are seen to give the same values when C ^ 3 x 10^"^ 
rad/s. 

The structure of Eq. (j3.ip shows that for numerical integration it is advantageous to 
introduce the nondimensional quantity 

y = qa (4.17) 
as the integration variable. The force expression then takes the form 



7T7, 6 ^ ^ B 6 ^ 



(4.18) 



(This formula holds even when practical units are restored, when j3 = l/fc^T.) Typical 
magnitudes of the attractive pressure are about one millipascal, for a gap width of 1 /xm. 
(The force between ideal metal plates at zero temperature for 1 /xm separation is 1.30 mPa.) 
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The next task is to determine the values of and Bm, in the hmiting case of m 0. 
This has to be done analytically. Whereas the TM mode leads unambiguously to Aq = 1 
{e ^ 1), the TE mode is more delicate. In Sec. Elwe show explicitly, by means of a limiting 
procedure based on the Drude model, how when ( ^ 0, i.e., when m — > 0. The 

m = TE mode accordingly does not contribute. To summarize: 

^0 = 1, Bo = for a metal (^(O) = oo), (4.19a) 

Aq = ^— j^Y^ ; Bq = for a dielectric medium {e = £{0)). (4.19b) 

These relations will be assumed in the following. 

There are some general properties of the expression ()4.18|) that ought to be noticed. First, 
at the lower limit, y = mj, the coefficients A^ and Bm for m > 1 become equal, 

Am = Bm= (^^^i^] , £ = £«m). (4.20) 



This expression is precisely the reflection coefficient for Poynting's vector, at normal inci- 
dence. This special case obviously corresponds to = 0. Then the TE and TM modes 
are identical to each other. Secondly, we note that for large values of y, the integrand in 
Eq. ()4.18|) approaches 

' e = e{^U, (4.21) 

showing how quickly the contributions from large y die out. 

The full line in Fig. |2] shows how the magnitude of JF^ for gold varies with the dimen- 
sionless parameter aT, when a = 1 /zm. The lower limit aT = 4.4 x 10~^ corresponds to the 
low temperature of T = 10 K. Terminating the y integration at the upper limit ?/max = 30 
we found the necessary number of terms in the m sum to be about = 450. At room 
temperature, T = 300 K, corresponding to aT = 0.131 and 7 = 0.823, the required number 
of terms was found to be lower, = 15 (assuming the same ?/max)- In the upper limit, 
aT = 0.52 (T = 1200 K), only = 4 was required. This property of only a small number of 
terms being necessary at high temperatures is as we would expect. Note, however, that the 
temperature variation of £{iC) is not taken into account. The only known empirical data for 
£{iC) are referring to room temperature, and are as given in Fig. ^ 

The broken line in the same flgure gives the result calculated from the expression in 
Eq. ()2.10|) . which is for the modifled ideal metal model in which the TE zero mode has 
been removed. The deviations from the full lines are seen to be quite uniform: 13% at 
the lower limit, 12% at room temperature, and 18% at the upper limit. This uniformity 
in the deviations is somewhat surprising, in view of the fact that the expression (j2.10p is 
a low-temperature expansion which one would expect to be most accurate when aT — 0. 
The reason for the deviations must lie in the different ways the two force expressions are 
calculated: Eq. ()2.10|) is based upon the idealized assumptions Am = Bm = 1 for all m 
except that Bq = 0, whereas Eq. ()4.18p is calculated using the reahstic dispersive data from 
Fig. 1, plus Eq. ()4.19ap in the case m = 0. 

Figure El shows that the behaviour is essentially the same if the gap is made wider, 
a = 4 fim. The forces are now only about 0.4% of those in Fig.|21 The lower limit aT = 0.017 
corresponds to T = 10 K (A^ = 115 terms necessary), and the upper limit aT = 0.523 
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a = 1 |im 




FIG. 2: Magnitude of surface force density for gold, in the temperature interval lOK <T < 1200 
K, when a = 1 ^m. Solid line is physical result calculated from Eq. (|4.18j) where the room- 
temperature data for e(iC) shown in Fig. ^ are used. Broken line is calculated from the ideal 
low-temperature form ()2.10|) . 

corresponds to T = 300 K (7 = 3.29, = 4). The deviations between the full dispersive 
result and Eq. ()2.10|) are now smaller than previously, about 5%. 

As experiments are usually made at room temperature for various gap widths, we show in 
Fig.EJhow the surface force density for gold varies with a, at T = 300 K. We have here chosen 
to multiply the ordinate with a^. The linear slope seen for a > 4/im is nearly that predicted 
in Eq. (j4.14p . which gives a slope of 2.0 x 10~^^ Nm^//im. The linear region between 1 and 
2 nm corresponds to that in Eq. ()2.10p or ()4.14|) (intermediate temperatures). Also shown 
is the prediction of the temperature dependent Drude model (Appendix^, when T = 300 
K. The differences are seen to be very small. Since the Drude values for the permittivity are 
lower than the empirical ones at high frequencies, as seen in Fig. Q we expect the predicted 
Drude forces to be slightly weaker than those based upon the empirical permittivities. This 
expectation is borne out in Fig. HJ the differences being large enough to be slightly visible at 
short distances, as we would expect since the plasma nature of the material becomes more 
pronounced for small distances. Note that the temperature dependence of the permittivity 
is irrelevant here because the temperature is fixed, unlike in Figs. |21 and IHl 

It is of interest to check the magnitude of the dispersive effect in these cases. We have 
therefore made a separate calculation of the expression ()4.18p when e is taken to be constant. 
Figure IHl shows how the force varies with aT in cases when e G {100, 1000, 10000, 00} are 
inserted in the expressions for Am and Bm in Eq. ()3.2a|) . Note that the e = 00 curve is 
obtained easily via the analytic result ()3.13|) . with = = 1 for all m > 1. With 
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a = 4 |am 



4.8 



4.6 



4.4 



(0 4.2 
Q. 

4 



I 



3.8 



3.6 



3.4 



3.2 



1 1 1 

\ 






\ 






\ 




- 


\ \ 






\ \ 






\ \ 






\ \ 






\ \ 






\ \ 






\ \ 






\ \ 






\ \ 






\ \ 




- 


\ \ 






\ \ 






\ \ 




/ 


\ \ 
\ \ 
\ \ 
\ \ 
\ \ 




/ 

/ 

/ 

/ 


\ \ 
\ \ 
\ \ 




/ 

/ 

/ 









0.1 



0.2 0.3 0.4 0.5 

aT 



0.6 



0.7 



FIG. 3: Same as Fig. |2l but at a larger spacing, o = 4 /im, corresponding to lOK < T < 300 K. 



T = 300 K 




4 5 6 

a(|am) 



FIG. 4: Surface force density for gold, multiplied with a^, versus a when T = 300 K. Input data 
for e(«C) taken from Fig. ^ 
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a = 1 |im 




FIG. 5: Nondispersive theory: Surface force density calculated from Eq. (|4.18|) for e G 
{100, 1000, 10000, cxd}. The £ = oo result is calculated from Eq. (|4.22|) . For low values of aT 
the latter coincides with the expression (|2.1())) used in Fig. |2I Also shown for comparison is the 
dispersive result for gold, where experimental input data for e(iC) taken from Fig. ^ Gap 
width is a = 1 /um. The constraint a = 1 fim applies only to the dispersive case, since otherwise 
a^T'^ is a function of aT only. 



i?o = 0, Eq. ()3.13p is modified into 

^^(e = oo) = ^ |c(3) - g 1 [s,{^k) + 2s,{^k) + 2.0(7^)] I , (4.22) 

which amounts to adding the last term of Eq. ()2.1()j) . [The sum is alternatively given in 
Eq. ()3.39|1 . and the low-temperature limit is given in Eq. ()2.1()|) .] It is seen from the figure 
that the three first curves asymptotically approach the e = oo curve, given by Eq. ()4.22|) 
when e increases, as we would expect. Again, we emphasize that the dispersive curve for 
gold is calculated using the available room-temperature data for e{iQ from Fig. [H In the 
nondispersive case, there is of course no permittivity temperature problem since e is taken 
to be the same for all T. 

There are several points worth noticing from Fig. El (i) The curves have a horizontal 
slope at T = 0. For finite e this property is clearly visible on the curves. This has to be 
so on physical grounds: If the force had a linear dependence on T for small T so would 
the free energy F, in contradiction with the requirement that the entropy S = —OF/dT 
has to go to zero as T ^ 0. For the gold data the initial horizontal slope is not resolvable 
on the scale of this graph, see the discussion at the end of Sec. IIVBI (ii) The curves show 
that the magnitude of the force diminishes with increasing T (for a fixed a), in a certain 
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y = l y = 3 



m e 


Xm) X 10^ Cm 


X lO^^rad/s 










1 


382.0 


8.226 


0.9998 


0.7899 


0.9999 


0.4944 


3 


100.4 


24.68 


0.9990 


0.8578 


0.9997 


0.6317 


5 


49.76 


41.13 


0.9975 


0.8774 


0.9992 


0.6759 


7 


30.28 


57.58 


0.9956 


0.8872 


0.9985 


0.6985 


9 


20.52 


74.03 


0.9931 


0.8930 


0.9977 


0.7124 


11 


14.87 


90.49 


0.9902 


0.8970 


0.9967 


0.7219 


13 


11.30 


106.9 


0.9867 


0.8998 


0.9955 


0.7288 


15 


8.891 


123.4 


0.9827 


0.9020 


0.9942 


0.7341 



TABLE I: Some data in the dispersive theory for gold. Here T = 10 K, y = go G {1,3}. Room 
temperature input data for e(«C) taken from Fig. ^ 



a(//m) 


m = 


m = 1 


m = 2 


m = 3 


m = 4 


m = 5 


m = 6 


m = 7 


0.5 


0.32 


0.98 


1.03 


1.05 


1.06 


1.07 


1.08 


1.08 


1 


0.58 


1.98 


2.05 


2.07 


2.08 


2.08 


2.07 


2.06 


2 


1.10 


4.04 


4.09 


4.07 


4.02 


3.96 


3.88 


3.79 


3 


1.63 


6.11 


6.09 


5.98 


5.80 


5.59 


5.36 


5.10 


4 


2.16 


8.18 


8.04 


7.75 


7.37 


6.93 


6.45 


5.95 


5 


2.69 


10.24 


9.92 


9.37 


8.69 


7.94 


7.16 


6.38 


6 


3.23 


12.30 


11.71 


10.81 


9.75 


8.63 


7.51 


6.45 


7 


3.78 


14.33 


13.39 


12.06 


10.55 


9.02 


7.56 


6.24 



TABLE II: Contribution from the various Matsubara frequencies for gold. What is given is the 
percentage of J-^ for each mode in the region m £ [0,7]. The temperature is T = 10 K. Room 
temperature input data for £{iC) are taken from Fig. ^ 

temperature interval up to aT ~ 0.3. This perhaps counterintuitive effect is thus clear from 
the nondispersive curves as well as from the dispersive curves in Figs. El and El (iii) It is 
seen that the curve for e = const. = 1000 gives a reasonably good approximation to the real 
dispersive curve for gold when a = 1 fim; the deviations are less than about 5% except for the 
lowest values of aT {aT < 0.1). This fact makes our neglect of the temperature dependence 
of e{i() appear physically reasonable; the various curves turn out to be rather insensitive 
with respect to variations in the input values of s{i(). (iv) One notes that the curves (for 
large e) in Fig. are consistent with the free energy ()4.14|1 using the rough approximation 
(14. 9p for Bm- Especially one notes the initial decrease of the magnitude of the Casimir force 
for increasing T when e is large. As discussed below Eq. ()4.12|) this is again connected with 
the counterintuitive negative contribution to the entropy, (v) Also, it can be remarked that 
i?o = is required when e is finite. Otherwise the curves in Fig. and thus the free energy, 
would have a finite slope at T = which again would imply a finite entropy contribution at 
T = in violation with the third law of thermodynamics. 

Instead of confining ourselves to a "black box" calculation of the force expression ()4.18p . 
it is desirable to break up the expression somewhat, to see how the various values of m 
contribute. We do this in Tables U Mllll for gold. The first two tables refer to the case T = 10 
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m = 


m = 1 


m = 2 


m = 3 


?Ti = 4m = 5m = 6m = 7 


0.5 


10.20 


31.24 


22.95 


15.09 


9.18 5.28 2.91 1.55 


1 


20.07 


49.37 


20.83 


6.97 


2.03 0.54 0.14 0.03 


2 


44.56 


49.87 


5.17 


0.37 


0.02 


Q 
O 


( u.yo 




U.OO 


U.Ul 




4 


88.88 


11.07 


0.05 






5 


96.58 


3.42 








6 


99.06 


0.94 








7 


99.76 


0.24 









TABLE III: Same as in Table Ull but at temperature 300 K. Data from Fig. ^ is again used. 

K. (Again, the experimental values of e{iQ at room temperature are used.) As y is the 
important integration parameter in Eq. ()4.18j) . we keep y fixed in Table HI y G {1,3}. It is 
seen that Am stays close to 1, whereas decreases for increasing y, if m is kept constant. 
Table HTl shows how the various m contribute to the force. Writing the total force as a sum, 

oo 

= E (4.23) 

m=0 

the columns in the table show the percentage of JF^, i.e., (JF^/JF^) x 100, distributed over 
the region m G [0, 7], when T = 10 K. The distribution from the various ms is seen to be 
very broad, as is characteristic for a low-temperature problem. Table IIIII shows the same 
kind of distribution over m when T = 300 K. Already from a gap distance of a = 3-4 fim 
onwards, the distribution is heavily concentrated around low m, as is characteristic of a 
high-temperature problem. 

It is in this context instructive as a corollary to go back to the integral over y in Eq. (|4.18|) . 
One would expext the main contribution to the integral to come from the region y = qa = 
^yk'j_ + ('^ a ~ 1. Assuming the most important values of k± to be moderate, this means 

~ 1, or m ~ l/(27raT), since ( = 27imT. When T = 300 K, we thus expect the dominant 
contribution to come from m ~ 1 when a = 1 fim, and from m = when a > 3 fim. This is 
seen to agree very well with the data in Table IIIII Similar considerations apply to the case 
T = 10 K, although the contributions from the various ms are then more smeared out. 

The important question is now: Have the characteristic temperature variations shown in 
the theoretical figures above been verified in practice? Of most interest in this context is 
the experiment of Bressi et al. [1^, since it deals with parallel plates directly. According 
to personal information from R. Onofrio, one of the members of the Italian group, the 
observed Casimir forces were lower than those predicted by the traditional (SDM) theory 
for conducting plates, in cases where the distances were low, a < 0.5 /im. This reduction 
effect is apparent also from their Fig. 4. Now, the plates in this experiment were coated with 
chromium rather than with gold, but we can check that the corrections in that case are of 
the same magnitude as if the plates were coated with gold. Namely, an explicit calculation 
of the analogue of Fig. El for the case a = 0.5 /zm (not shown here) shows that at room 
temperature for which aT = 0.065, the force becomes —J-^ = 15.5 mPa. The conventional 
(SDM) theory gives in this case the force 1.3 x 2^ = 20.8 mPa. The predicted reduction in 
the force is thus about 25%, somewhat more than the measurements indicate. In any case. 
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this suggests that the reduced force seen at room temperature in Ref. j26| may be the first 
actual observation of the temperature effect predicted theoretically. 

At larger distances, however, between 1 and 2 /xm, the situation is no longer so clear-cut, 
since they observe a Casimir force in excess of the theoretically predicted one. The reason 
for this deviation is not known. Of course the force becomes weaker at larger distances, 
thus being subject to larger experimental uncertainties. The most natural conclusion to 
be drawn at this stage is that we have to wait for better precision in this kind of difficult 
experiment. Ideas for such an improved experiment which could descriminate between the 
different models have just appeared ji^. 



APPENDIX A: ON THE SMOOTHNESS OF THE REFLECTION COEFFI- 
CIENT r2 AT LOW FREQUENCIES, FOR A METAL 

In view of the current discussion in the literature about the value of the reflection co- 
efficient r2 for a metal in the limit of low frequencies, let us consider this point in some 
detail. As mentioned earlier, the problem occurs in connection with use of the Drude for- 
mula, Eq. (jl.lj) . The coefficient r2 is actually the square root of our quantity defined in 
Eq. ()3.2a|) . so that we may write 

/ \ 2 

2 I S-P 



Let us keep the transverse wave vector fixed, and perform a power series expansion of 
s{iC,) to the first order in C/z/. [Any normal metal must have a finite relaxation frequency 
so that in the limit of low frequencies, C,/v can be regarded as small. At zero temperature, 
we are assuming v{T = 0) 7^ 0.] From Eq. (jl.ip we get 

e(zC) - 1 - ^ fl - -V (A2) 



which for the Lifshitz variables s and p implies [cf. Eq. ()3.2b|) ] 



c 

Insertion into Eq. (jAlll now yields 



,.^^3TT7^^(l + ||), (A3a) 
P = ^,/l + S-¥- (A3b) 



{-] ■ (A4) 



2 
P 

AklJ \u 

We thus see that — smoothly as C — > 0. Contrary to recent statements in the literature 
wc find that there is no peculiar effect taking place at C = 0, when the Drude model 
is used. The result flA4jl corresponds to a vanishing contribution to the Casimir effect from 
the m = TE mode for a real metal, in accordance with our treatment in Sec. IIVI 

The argument above hinged on the assumption that 7^ 0. One might wonder: What 
happens if is exactly zero? Mathematically, it then follows from Eq. ()3.2b|) that r| = 1. 
This case cannot, however, be of physical importance. The set k^ = is mathematically of 
measure zero, and has thus no influence upon real physics. 
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APPENDIX B: ON THE PHYSICAL IMPORTANCE OF AND 

It is physically instructive to show in some detail how the coefficients Am and Bm relate 
to the conventional Fresnel coefficients in optics, at oblique incidence. Consider ffist the TM 
mode, and let a plane wave be incident from the left (medium 1, refractive index rii = ^/6) at 
a real angle of incidence 6i towards the boundary located at z = 0. The angle of transmission 
to the vacuum region z > is 6t. For instance from Ref. ||47| we have for the ratio between 
the reflected wave amplitude and the incident wave amplitude 

i?™ cos — 77-1 cos , . 

A™ ^ cos^, + ni cos^t' ^ ^' 

Since cos = a/1 — k^/ (£<^^), cosOt = a/1 — k^/uj"^ we get, when replacing uhy iC,, 

R™ _ ^e + ki/e-E^i + ki/e 



^™ V^ + kl/C + e^l + kl/C' 



(Bib) 



Now s = ^/ e — 1 + p'^ = a/^ + kj_/C,'^, P = q/C = a/1 + k^/C^, and so we get 

R™ s-ep 



A™ s + ep 

Similarly for the TE mode. 



'Am. (B2a) 



^™ ' ^ (B2b) 



Of course, these results are also found in textbooks |4 
APPENDIX C: PARALLEL DIELECTRICS 

In Ref. P the following result for the TE reduced Green's function is given, 

g"{z, z') = — f e-"2l"-"'l + r e-''2{.+.'-2a)^ _ ^^^^ 



2K2 

which is valid for z, z' > a. Here the reflection coefficient is 



r = + ^d-\ (C2a) 

d = ^l±^^l±^e''^^^-l, (C2b) 

f^3 - f^lK'3- f^2 

and 

^2 = k'^ - ^^.^ (C3) 
and we have taken a parallel dielectric slab geometry 

ei, z <0, 

e(^) = { e-3, 0<z<a, (C4) 
€2, a < z. 
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The temperature controversy centers on the zero mode. If uJ^e vanishes at = (true for 
the Drude model, but not the plasma model), then the reflection coefficient vanishes there, 
r = 0, and we have only a free Green's function at a; = 0, that is, the boundary becomes 
transparent. The TM reflection coefficient does not have this property. 

We have redone the calculation to find the reduced Green's function in the interior region, 
^ < z,z' < a. We find 



g'^iz, z') 



1 

2^3 I 



K3 — Ki 

( 

K3 + Ki 



-K:i[z-^z') _|_ ^-1 



gK3(2-z') _|_ ^Kz{z'~z) 



K3 - Ki fi;3 + Ki 

Again, it is easy to see that we obtain only the free Green's function for the zero mode 

1 



(C5) 



g^{z,z']u; 



0) 



-k\z~z'\ 



2k 



(C6) 



provided lim^^o = 0. 

A check of this result is that if we substitute Eq. ()C5|) into the expression for the force/area 
(3.13) of Ref. yj, we get for the TE contribution to the force [see Eq. (3.10) there] 



2ti (27r)3 
duj ((ik) 
2^(27r; 



62 - e3)uj^g^{a,a) 
3 (/t3 - ^2 + 2Ksd~^) 



(C7) 



identical to the first term in Eq. (3.19) of Ref. [l|, and apart from a contact term is the 
same as the second t^erm in Eq. ()3.1|) . All of this does not seem to support the claims of 

|38| that there is something ill-defined about the u; = limit. 
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APPENDIX D: TEMPERATURE DEPENDENCE OF THE RELAXATION FRE- 
QUENCY FOR GOLD 

To investigate the temperature dependence of the relaxation frequency z^(T) in the Drude 
relation 

for gold, it is convenient to make use of the Bloch-Griineisen formula for the temperature 
dependence of the electrical resistivity p ji^ : 

It is known that = 175 K for gold. The constant C can be determined from the knowledge 
that p = 2.20 X 10-^(] m at temperature 295 K [501 . We obtain C = 5.32 x lO^^r^m. 
The theoretical relationship between i> and the static resistivity p is 

u = p, D3 

m 
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300 



T(K) 



FIG. 6: Temperature dependence of the relaxation frequency for gold. 



where N^. is the number density of atoms, foN^. with /q ~ 1 the number density of free 
electrons, and m the effective electron mass. The simplest way to proceed is to put v = Kp 
with K a constant, and make use of the room-temperature data of Eq. ()4.16|) . We obtain 
K = 1.59 X 10^ eV n-^m-\ Altogether, 



z/(T) = 0.0847 



(e- - 1)2 



(D4) 



where the unit of i^(T) is eV. The temperature variation is shown in Fig. IHl For low tem- 
peratures, i^{T) oc T^, whereas at high temperatures, i^{T) oc T. The curve is seen to be 
similar to the one given in Fig. 3 of Ref . , in the case of aluminum. 

An important caveat must be mentioned, however; these formulas ne glec t the effect of 
impurities, which give rise to a nonzero resistivity at zero temperature j5l|. This makes 
the use of these ideal resistivity models questionable, and adds further evidence that the 
behavior of the entropy discussed in Sec. |TV|is correct. 
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